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ABSTRACT 

We use a combination of reverberation mapping data and single-epoch spectra of the C IV emission 
line in a sample of both low and high-redshift active galactic nuclei ( AGNs) to investigate sources of the 
discrepancies between C iv- and H/3-based single-epoch black hole mass estimates. We find that for all 
reverberation mapped sources, there is a component of the line profile that does not reverberate, and 
the velocity characteristics of this component vary from object-to-object. The differing strength and 
properties of this non- variable component are responsible for much of the scatter in C iv-based black 
hole masses compared to H/3 masses. The C iv mass bias introduced by this non-variable component 
is correlated with the shape of the C iv line, allowing us to make an empirical correction to the black 
hole mass estimates. Using this correction and accounting for other sources of scatter such as poor 
data quality and data inhomogeneity reduces the scatter between the C IV and H/3 masses in our 
sample by a factor of ~2, to only ~0.2 dex. We discuss the possibility that this non- variable ClV 
component originates in an orientation-dependent outflow from either the proposed broad line region 
(BLR) disk- wind or the intermediate line region (ILR), a high- velocity inner extension of the narrow 
line region (NLR). 

Subject headings: galaxies: active — galaxies: nuclei — line: profiles — quasars: emission lines 



1. INTRODUCTION 

Observations of galaxies in the local universe have led 
to the nearly universal acceptance that massive galax- 
ies house a supermassive black hole (BH) at the cen- 
ter of their gravitational potential well. In the local 
universe, the masses of the BHs are strongly corre- 
lated with some properties of their host galaxies, lead- 
ing t o the Mrh — o* and Mrh — -khnire relationships 
(e.g. IFerrarese fc Merrittl l2000t iGebhardt et~al1 120001: 
Tremaine et al.l I2002t Marconi fc Huntl 120031 ; IGrahaml 
2007tlGultekin et al.ll2009HGraham et al.ll2011h . Because 
of the apparent tight evolutionary relationship between 
the BH and its host, one way in which to probe galaxy 
evolution is to trace the growth of the central BH. In 
order to do this, we need reliable methods for measuring 
BH masses in the near and distant universe. 

Direct dynamical BH mass measurement methods 
work only for nearby quiescent galaxies in which the 
gravitational sphere of influence of th e BH can be spa- 
tially resolved (see, e. g., the review by IFerrarese fe ForcH 
120051 though also see IGiiltekin et al.l 120091 for a discus- 
sion of deriving masses without resolving the sphere of 
influence). On the other hand, reverberat ion mapping 
(jBlandford fc McKed [1981 iPetersonl [1991 can be ap- 
plied to active galactic nuclei (AGNs) to measure BH 
masses both locally and at cosmological distances. Un- 
fortunately, the required length of the spectrophotomet- 
ric monitoring campaigns needed to make the measure- 
ments becomes longer with increasing source luminosity 
and distance. As a result, current mea surements only 
extend to sources with z < 0.3 (see, e.s? iPeterson et al 



2004; 



Bentz et a l. 2009bl; iDennev et al.feoiOtlGrier et a! 



2012at see Kaspi et al. 2007 for a tentative result for a 
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z ^2 quasar). 

Extending direct mass measurements to large samples 
of objects at high redshifts in order to study the co- 
cvolution of BHs and galaxies is not likely to be pos- 
sible in the near future with current resources. However, 
the results from reverberation mapping of local AGNs 
can be used to calibrate indirect estimates of BH masses 
for large samples of objects using single spectra. These 
"single-epoch" (SE) mass estimates rely on the observed 
tight correlation between the monochromatic AGN lu- 
minosity, L, and the radius of the broad line region 
(BLR), R, measured in reverberation mapping experi- 
ments. With this method, only two observables are then 
needed: (1) a measurement of the AGN luminosity to 
use as a pro xy for the BLR radius through this R — L 
relatio nship (jKaspi et al.l 120051 120071: iBentz et al.l 120061 
2009a), and (2) a measurement of a broad emission- line 
width (typically Ha, H/3, MgnA1549, or CivA1549) to 
estimate the Doppler-broadened BLR line-of-sight (LOS) 
velocity dispersion, V. Combining these observables, a 
virial BH mass can be calculated as Mbh = fRV 2 /G, 
where G is the gravitational constant and / is a scale 
factor of order unity related to the BLR geometry and 
kinematicfl 

The majority of reverberation mapping results to 
date use time delays measured from H/3 emis- 
sion line variability. Consequently, the most ro- 
bust SE mass scaling relation is based directly 
on the R — L relat i onship for this emission-line 
(see ICollin et al.l 120061: IVestergaard fc Peterson! 120061: 



3 For reverberation mapping results that use the line width 
measured from the rms spectrum, an ensemble average, {/), sets 
the AGN BH mass scale zero-point, and is determined assuming 
AGNs follow the same M — a* rela tionship as quiescent galaxies 
HOnken et al.ll200l IWoo et al.ll20TOl) . For SE masses, / has a dif- 
ferent calibration based on the differences between the SE and rms 
line widths for H/3 (see lCollin et~aT|[200rD . 



2 



iBentz et al.l [2009a). Scaling relatio nships for other 
broad emission lines are also availab l e dMcLure fe Jarvisl 
2001 IVestereaard fe Petersoiil [20061: IMcGill et al.l 120081: 
Vestergaard fe Osmeri 120091 ) . but these SE mass scaling 
relationships are calibrated to the RM-based H/3 mass 
scale rather than by direct measurements of the time 
delays for these lines. In this way, SE masses can be 
estimated for large samples of AGNs observed in spec- 
troscopic surveys and across a large range of redshift s 
(see, e.g- IKollmeier et al.l l2006: Shc n et alJl2(M [201ll) . 

This suite of intercalibrated AGN BH mass scaling re- 
lations initially seems the ideal method for using the 
myriad of optical survey spectra of AGNs to study BH- 
galaxy coevolution at all redshifts. However, there are 
significant concerns as to the reliability of masses based 
on lines other than H/3 due to observed systematics 
in these masses compared to H/3 masses. For exam- 
ple, only tentative re verberation results exist for Mgn 
(|Metzroth et al.lfeOOGD . so the use of Mg II as a virial mass 
indicator is based only on its simila r ionization poten- 
tial and line width compared to H /3 (jMcLure fc Dunloq 
12004 IVestereaard fc Osmeri l2009f ). Studies have shown, 
however, that Mgn-based masses suffer from clear sys- 
temat ics, some of which may be related to the Eddington 
ratio (Onk en fc Ko llmcicr 2008) or syste matic line widt h 
differences between Mgn and H/3 (see iCrooml feOllf) . 
Mgn-based masses are also susceptible to biases related 
to the prescription used for measuring the line width 
(jRafiee fc Halil20Tl . 

The same is true of SE masses based on 
C IV. The current C IV m ass scaling relationship 
(jVestergaard fc Peterson|[2006l hereafter VP06) is based 
on a direct calibration to H/3 reverberation-mapped 
AGNs. For this sample, VP06 found a relatively 
small scatter in the SE masses compared with the 
RM masses, ~ 0.3 dex, which has since been cited as 
the typical assumption of the accuracy of SE masses. 
Additional studies have also found a general consistency 
between Civ-based and B aimer- line- b ased masses (e.g., 
iGreene et all 120101 : lAssef et al.1 120111 ). However, other 
authors have questioned the reliability of C iv as a result 
of finding little consistency or correlation and a large 
scatter between C iv- and H/3-based SE masses. These 
studies contend that Civ-based masses have too much 
scatter compared wit h H/3 masses to be a reliable virial 
mass indicator (e.g.. iBaskin fc Laorll2005|: [Netzer et al.l 
l200l [Suientic et all 120071 : iShen fc Liull2012l ). 

One common postulate for this seeming unreliability 
is that the C iv em ission region is susc eptible to out- 
flows and winds (see lRichards et aLll2011l and references 
therein) . Such non- virial gas velocities are then believed 
to bias the resulting SE mass measurements, rendering 
them unreliable. Interestingly, however, such dynam- 
ics do not seem to affect the C iv RM-based masses. 
In the few objects where RM measurements have been 
made for multiple emission lines, C iv results follow the 
expected virial relation w ith other low-ionization lines 
(|Peterson fc Wandell [20001 and the BH masses derived 
from the i ndividual lines, inclu ding C IV, are mutually 
consistent ([Peterson et al.l [2004) . Here, we explore this 
apparent contradiction — why would RM-based Civ 
masses behave as expected but SE masses do not? In $2] 
we first review the fundamental assumptions that make 
estimating a SE mass possible. In Sj3]we investigate prop- 



erties of the C iv profile in the RM sample, identifying 
a source of bias related to the Civ line profile. In 2] 
we quantify this bias and empirically fit a correction to 
C IV masses that significantly reduces the scatter between 
Civ and H/3 SE masses. Finally, in §$5]and[6]we put our 
results in context to other studies, suggest a physical in- 
terpretation to explain our observations in terms of an 
orientation-dependent outflow, and summarize our find- 
ings. 

2. SINGLE-EPOCH MASSES: THE FUNDAMENTAL 
ASSUMPTIONS 

We start by reviewing the assumptions made to esti- 
mate a SE BH mass. These assumptions are based on 
the different means with which the physical parameters 
required to measure a black hole mass (BLR radius and 
velocity) are inferred in indirect (SE) and direct (RM) 
methods. 

SE masses rely on the continuum luminosity as a proxy 
for the BLR radius, and this is susceptible to two sources 
of uncertainty that do not affect direct RM measure- 
ments of the BLR radius: variability and host galaxy 
contamination. First, since RM measures the BLR ra- 
dius directly from the delay between the continuum emis- 
sion and the reprocessed line emission, an RM mass is not 
directly dependent on the luminosity. Instead, the mea- 
sured BLR radius and velocity (and associated uncertain- 
ties) are physically connected to the intrinsic AGN lumi- 
nosity through photoionization and its effects on BLR 
emission pro perties (e.g., the loca lly optimally-emitting 
cloud modeh lBaldwin et aDll995f ). As a result, RM will 
always give th e same mass independent of the intrin sic 
variability (cf. IBentz et all 120071 : iGrier et al.l l2012allbl ). 
However, the delay that makes RM experiments possible 
means that the photons responsible for the luminosity 
measured in a SE spectrum are not the same as those re- 
sponsible for producing the line emission whose width is 
measured in the same spectrum. Therefore, the reliability 
of SE masses depends on the assumption that the differ- 
ence between the observed SE luminosity and that respon- 
sible for triggering the observed line emission is negligi- 
ble, or at least small. Fortunately, the variability ampli- 
tude of AGNs is relatively small, ~ 10—30%, on reverber- 
ation time scales — the li ght-cr ossing time of the BLR 
(see, e.g.. iMacLeod et al.l f20ToT ) . Since M SE oc L~ - 5 , 
the uncertainty in the luminosity due to variability is 
not a dominant source of uncertainty in SE masses and 
is largely accounted for in the error budget through the 
uncertain ties in the calibrati on of the R — L relationship 
itself (see IBentz et al.ll2009af) . Second, only the AGN lu- 
minosity is correlated with the radius of the BLR, since 
photoionization of the BLR gas from the AGN continuum 
photons is what drives the R — L relationship. There- 
fore, SE masses will be overestimated if the observed 
SE luminosity is not corrected for host galaxy starlight 
contamination. This contamination has been removed 
from RM data sets used to calib rate the R — L rela- 
tionship ([Bentz et al.ll2006l [2009a|) . but the RM masses 
themselves are independent of this contamination. 

Next, because velocity is as important a p arameter as 
the luminosity in determining SE BH masses ( Ass ef et al.l 
|20T1 and M oc V 2 , we must consider differences in the 
BLR velocity field characterizations of the SE and RM 
methods. RM masses are largely based on line widths 



CIV Line Shape 



3 



measured from the rms spectrum (i.e., the spectrum 
formed from the rms flux deviations about the mean 
spectrum, based on all spectra obtained during a sin- 
gle reverberation mapping campaign), so that only the 
kinematics of the gas reverberating in response to the 
continuum photons contribute to the mass determina- 
tion. The SE line profile is not necessarily equivalent to 
the rms profile because there can be non-variable emis- 
sion components present. So, the reliability of SE masses 
also depends on the fundamental assumption that the SE 
line-width is an accurate proxy for the rms BLR veloc- 
ity field. Comparisons of mean and rms H/3 emission-line 
profiles show that even after removing the contamination 
from narrow line region (NLR) emission, the SE pro- 
file is not identical to the rms profile — not all the H/3 
broad line-emitting gas is reverberating, particularly in 
the high- velocity wings. This results, on average, in nar- 
rower rms spectrum profiles than mean spectrum profiles. 
For H/3, these differences are not significant enough to 
strongly bias the SE mass estimates once the differences 
are taken into account statistically by the c alibration of 
the sc ale factor, ( f), fo r SE masses (see, e.g.. lCollin et al.l 
l2f)M IPark el~allf20ll . 

3. THE SE VERSUS RMS C IV PROFILE 

The many studies comparing SE H/3 masses to RM 
masses reinforces the reliability of these masses, but no 
similar investigation exists for the other emission lines 
used for SE masses. The implication of (1) the large 
scatter between C iv and H/3 SE masses and (2) the con- 
sistency of C iv and H/3 RM masses strongly suggests 
that the problem lies in the assumption that the C iv 
mean and rms line widths are similar. 

To examine this, we collected mean and rms line width 
characterizations from RM experiments of both H/3 and 
Civ. We take the H/3 FWHM and line dispersion, 
o~i (the second momen t of the line profi le) , measure- 
ments from Table 1 of iCollin et al.l ()2006l ) . which have 
all been me asured homogeneously using the methods de- 
scribed by iPeterson et al.l (|2004| ). Such a resource is 
not available for the C IV RM data sets, so we acquired 
the original mean an d rms spectra from all available 
Civ RM experiments dClavel et alj|1991t iReichert et al.l 
1994HKorista et al.lll995HRodriguez-Pascual et al.lll997l; 
Wanders et al.lll997t lO'Brien et al.lll998t IPeterson et al l 
20051: iMetzroth et al.ll2006t iKaspi et al.l l2007f T From 
these mean and rms spectra, we fit and subtracted a 
linearly interpolated continuum between restframe con- 
tinuum regions near 1450A (except for the NGC 7469 
rms spectrum, which was fit near 1500 A, where the con- 
tinuum appeared lower) and 1700A. We then measured 
the Civ line widths directly (i.e., no functional forms 
were fit), except for the case of the rms spectrum of NGC 
4395, which was too noisy to measure directly. In this 
sixth order Gauss-Hermite polynomial was fit to 
the data, and the width was measured from the fit. We 
employ the same line width measurement techniques as 
those utilized for H/3 in order to have a homogeneous 
comparison. These new C iv measurements are given in 
Table 1. 

Figure [T] shows the ratios of the mean and rms line 
widths as a function of the line width for H/3 (open 
symbols) and ClV (closed symbols). We consider first 
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Figure 1. C IV and H/3 line width ratios measured from the mean 
and rms spectra from reverberation mapping experiments. The 
left and right panels show ratios of the widths measured using the 
FWHM and ct;, respectively. Solid (open) points represent mea- 
surements of Civ (H/3) emission-line widths. The solid horizontal 
lines corresponds to a ratio of unity. 

the results for the FWHM. The known trend that H/3 
shows systematical ly narrower rms than mean line widths 
(|Collin et al.l [2006) is clearly apparent, as the H/3 ra- 
tios lie mostly below unity. On the other hand, C iv 
shows a systematically different trend, with all objects 
having broader rms than mean spectrum line widths, 
and in some cases the difference is quite large. The 
line dispersion ratios for H/3 show a similar offset from 
unity, albeit with greater scatter. On the other hand, 
the C iv line dispersion ratios are now uniformly scat- 
tered around unity. Figure Q] shows that there are signif- 
icant differences between the time variable and instanta- 
neous or mean characterization of the H/3 and C iv broad 
emission-line widths. These differences are larger for the 
FWHM than for the line dispersion, and they are the 
largest and least uniform for the FWHM characteriza- 
tion of the C iv line. 

We investigate this incompatibility further by parame- 
terizing the shape of the line profile in terms of the ratio 
S = FWHM/cr;. Profiles with lower values of S have 
'peaky' profiles with larger kurtosis, while those with 
higher values of S have 'boxy' profiles with smaller kur- 
tosis. A Gaussian profile has S = 2.35 and a kurtosis of 
three. Figure [2] shows the ratio between the line shapes 
determined from the mean and rms spectra, S rms / S mean , 
for both H/3 and C iv. The shapes of the mean and rms 
H/3 line profiles are essentially indistinguishable. There 
may be a weak effect associated with narrow line Seyfert 
l's (open triangles), where it is difficult to accurately re- 
move the contribution from the narrow-line component. 
For Civ, however, the shapes of the mean and rms line 
profiles tend to be very different. The consistency in the 
mean and rms H/3 line shape allows the simple offset in 
H/3 widths seen in Figure Q] to be easily compensated for 
in the standard recalibration of the RM BH masses for 
SE data. No such calibration would work for C iv with- 
out introducing extra scatter in the SE mass estimates 
compared to RM mass estimates. 

The differences between the SE and rms C iv profile 
can also be seen through visual inspection of the mean 
and rms spectra observed in C IV reverberation mapping 
campaigns. Figure [3] shows the mean and rms profiles for 
th e RM sample, ex cluding the z ~ 2 quasar monitored 
bv lKaspi et ail (pOOl . The black curves show the mean 
spectrum C iv profile, normalized to the peak flux after 
subtracting a local, linearly interpolated continuum, and 
the gray, dashed-dotted curves show the corresponding 
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Table 1 

C IV Emission Line Widths of Reverberation Sample 



mean spectrum 



rms spectrum 



Object 


Data Set 1 


FWHM 




FWHM 




Ref.^ 


3c390.3 




6000±150 


5360±40 


11340±700 


5760± 50 


1 


Fairall 9 




3040± 50 


4100±40 


3290±900 


2670±130 


2 


NGC 3783 




3120± 60 


2910±20 


4060±200 


2790± 70 


3 


NGC4151 


1988 


3180±130 


5190±30 


13820±530 


5740± 70 


1 


NGC 4151 


1991 


6610±260 


4800±30 


10450±280 


5160± 90 


1 


NGC 4395 


Visit 2 


1070± 10 


2130±20 


5840±350 


2510± 50 


5 


NGC 4395 


Visit 3 


1270± 10 


1950±20 


1900±170 


990±150 


5 


NGC 5548 


HST 


4080± 70 


3940±10 


6840±150 


3640± 30 


6 


NGC 5548 


IUE 


4860±250 


4100±60 


5590±200 


4110± 60 


7 


NGC 7469 




3290± 90 


3570±30 


4380±190 


2670± 30 


8 


S50836+071 3 




8920±190 


4600±30 
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1 For objects with multiple measurements, we include 
for details. 

2 Orig inal Data References: (IjJO'Bricn ct al. (1998); 



(4) IMetzroth et al.) l20Qi 
(8) IWanders et al.l U997) 



(9) Kaspi ct al. (200 



note to distinguish between them; see References 

fgl lRodriguez-Pascual et al.l 1T99 1 
0] (|207~ 



(3) Rcichcrt et al. 
|); (7) IClavel et al.1 



3 Unfortunately, we had to exclude this object from our analysis because the rms spectrum was very noisy, 
and the Civ variable emission was blended with variable Hell A 1640, so a reliable Civ line width could not 
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Figure 2. Civ and H/3 line shape ratios measured from the mean 
and rms spectra from reverberation mapping experiments. Solid 
points represent measurements of C IV emission-line widths, open 
points refer to H/3 widths. The H/3 measurements have been di- 
vided by symbol type into those from normal broad-line AGNs 
(open squares) and NLSls (open triangles). The solid horizon- 
tal line represents the 1:1 relation of equal mean and rms shapes. 
(There are no Civ points with SWan > 1.4). 

continuum-subtracted normalized C IV line profile from 
the rms spectrum. Here the continuum was fit to both 
the mean and rms spectra using regions near restframc 
1450A and 1700A. Normalizing in this way makes it easy 
to see that the FWHM in the rms spectrum is broader 
than that in the mean spectrum for all objects, as was 
also seen in Figure [T] More interestingly, however, are 
the visibly different shapes of the C IV lines in the mean 
and rms spectra. This was also already demonstrated in 
Figure [5J but we see in Figure [3] that the mean spectrum 
has an additional, low-velocity component absent from 
the rms spectrum for all objects but Fairall 9. In the case 
of Fairall 9, the rms signal is weak and the strength of 
the narrow core component in the rms spectrum is within 



the level of the flux uncertainties of the IUE campaign 
data, which makes it difficult to interpret t he shape dif- 
ferences (see iRodriguez-Pascual et al.lll997[ ). In general, 
it appears that the SE or mean C iv profile is made up of 
two components: a non- variable, largely core, component 
plus a variable component, the velocity profile of which 
can be isolated in the rms spectrum. We can empha- 
size which portion of the mean spectrum is variable with 
the solid gray curves in Figure [3] that show the normal- 
ized rms spectra scaled in flux by an arbitrary fraction 
to match (by eye) the mean and rms spectrum flux level 
in the red wing of C IV. 

4. SHAPE BIAS IN C IV BH MASSES 

The existence of two components in the SE C IV pro- 
files, only one of which is variable, has significant implica- 
tions for the determinations of SE C IV BH masses, both 
for individual estimates but also for calibrating the C IV 
mass scale itself. The contribution of the non-variable 
component creates a bias to the relevant (i.e., reverber- 
ating) BLR velocity field. This is particularly true when 
using the FWHM to characterize the line width. Fur- 
thermore, the strength and profile of the non-variable 
component is not uniform object-to-object, so a simple 
change in the scale factor, (/), or subtraction of a con- 
stant, narrow velocity component will not work. Instead, 
it seems likely that any correction must be based on the 
shape of the individual C IV line profile. 

To explore this, we compiled several samples from the 
literature that (1) have both C IV and H/3 mass estimates, 
and (2) for which both FWHM and ai characterizations 
of the line width are published or for which we can mea- 
sure them from the original C iv spectra so that we can 
calculate the shape parameter, S = FWHM /a;. These 
samples include those of VP 6. lAssef et al.1 (|2011l here- 
after All), and lNetzer et al.1 (|2007l hereafter N07). All 
provide all the needed quantities in their Tables 3 and 4. 
We use their "Prescription A" because it is most similar 
to the method used by VP06. We also consider both the 
group I and group II objects from All, excluding only 
the broad absorption line quasar, H1413+117, leaving a 
sample of 11 objects. 
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Figure 3. Mean (black) and rms (gray dot-dashed) spectra of the 
C IV reverberation mapping sample. Spectra have been continuum 
subtracted and normalized in flux based on the Civ emission-line 
peak. Solid gray curves show the rms spectra scaled by an arbitrary 
factor to approximately match the red wing flux between each rms 
and mean spectrum (see Section \3\ for details). The red curves in 
the panel for NGC 4395 show the normalized (red dot-dashed) and 
scaled (red solid) Gauss-Hermite polynomial fit to the rms spec- 
trum used for line width measurements. We show only one set of 
spectra for each object, even though some objects were monitored 
more than once; results not shown are qualitatively similar to those 
here. 
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Figure 3. Continued. 

The VP06 sample is of particular interest for investi- 
gating biases, since it calibrates the C iv SE BH mass 
scale. VP06 often used multiple SE spectra of a single 
object and calculated masses from each epoch. Further- 
more, the epochs and number of spectra they used varied 
from object to object and between C IV and H/3. Here we 
are more interested in the SE C iv and H/3 masses and S 
for different objects rather than different epochs, so we 
derived a geometric average of each of these quantities 
for the VP06 sample. For objects with more than one 
epoch of data, we use the C iv FWHM and 07 to calcu- 
late log(S) for each epoch, and then calculate a simple 
average of the individual values of log(S'), adding uncer- 
tainties from each epoch in quadraturfl We also aver- 

4 Intrinsic variations in 5" could occur from epoch to epoch, so 



age the SE C IV masses and uncertainties given in Table 
2 of VP06 for each object in the same way. Finally, we 
recalculate the SE H/3 mass using the updated formula 
provided by All (their Equation 4), which is based on 
the most recent cal ibration of the H/3 R — L relation of 
iBentz et~a l. (2009a). We then similarly average these SE 
masses to end up with one H/3 mass for each object. Ob- 
jects that had SE H/3 masses but not SE C iv masses were 
omitted from the sample. For the six objects lacking an 
SE H/3 mass, we used the RM mass, which is generally 
equivalent^. This results in a sample of 27 objects taken 
from VP06. 

N07 did not publish er; values for their sample, but their 
study was based on SDSS spectra. Thus, we obtained the 
SDSS spectra and remeasured both the FWHM, 07, and 
their uncertainties following the methods for "Prescrip- 
tion A" described by All. This includes fitting sixth 
order Gauss- Hcrmite polynomials to the C iv profiles be- 
cause of the typically poor data quality of this sample 
(average S/N per pixel in the continuum near restframe 
1700A is only 5.4). Interestingly, in several cases, our 
measurement for the FWHM of the C iv line in the N07 
sample differs significantly from the values presented by 
N07, as shown in Figure 2J Unfortunately, N07 neither 
describe how they determined the widths of the lines 
given in their Table 2 nor provide any uncertainties on 
these measured values, so reconciling these differences is 
not possible here. Nonetheless, this highlights the dif- 
ficulties in these types of studies as different measure- 
ment techniques applied to the same data set can re- 
sult in statistically different measurements of the same 
quan t ities (see also lAssef et al.l 120111 : iVestergaard et all 
1201 U IPark et all 120 111 ). We will continue to use our 
own FWHM measurements of the N07 sample because 
they have been determined with methods consistent with 
those applied to the other samples used here. We flag 
seven of these objects not necessarily because they are 
all outliers in any part of our analysis, but because we 
have evaluated their line widths to be potentially un- 
reliable due to evidence for absorption in the Civ line 
profiles of the SDSS spectrfl 

Using our C IV line widths and the continuum lumi- 
nosities given by N07, we calculate new SE C iv masses 
for each of the 15 objects in this sample based on the 
FWHM and 07 following Equations 7 and 8 of VP06 (also 
Equation 6 of All). We also recalculate the SE H/3 mass 
using the FWHM and luminosity provided by N07 but 
using the formula given by All. Unfortunately, we could 
not calculate uncertainties in the mass estimates because 
N07 do not provide uncertainties on the H/3 width or 

averaging over epochs may be removing information about long 
timescale changes in emission-line profiles known to exist. How- 
ever, the standard deviation about the mean S in all objects for 
which we calculate an average is much smaller than the formal 
uncertainties in S propagated from the uncertainties in the mea- 
surements of the SE FWHM and <t; . 

5 This is not strictly correct because host galaxy starlight has not 
been removed from the SE luminosities. The resulting H/3 masses 
for these six objects will be somewhat underestimated compared 
to the expected SE H/3 mass (this underestimation is likely < 0.2 
dex, based on the average difference calculated from the rest of the 
sample, for which both SE and RM masses are available). 

6 The classification of absorption here is based solely on the 
observed SDSS spectra. However, the existence of some such ab- 
sorption has been verifie d in new, higher S/N data obtained for 
seven of the N07 sources ( Assef et al. 2011, Denney et al. in prep). 
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Figure 4. Comparison of Civ FWHM measurements taken from 
N07 for their sample of 15 SDSS AGNs and those remeasured in 
this study. Error bars are only shown for our new measurements 
because N07 do not provide uncertainties on their width measure- 
ments. The solid line shows the expected 1:1 relation in width. 

optical/UV luminosity measurements. Nonetheless, this 
partial reanalysis places all masses on the same mass 
scale, with homogeneously measured C IV line widths. 

In Figure [SJ we show the ratio of C iv to H/3 SE masses 
as a function of the line shape, S. The top panel shows 
the BH masses calculated using the FWHM of the line. 
The different symbol types represent the different sam- 
ples. There is a statistically significant correlation of the 
mass residuals with C iv line shape. The Spearman rank 
order coefficient, given in the top right corner of the plot, 
is r s = 0.601, and the probability that such a correlation 
would arise by chance if the quantities were intrinsically 
uncorrelated is P ran = 9.93 x 10~ 6 . These statistics were 
determined with only 46 of the 53 original objects be- 
cause we omitted the seven absorbed N07 targets (gray 
stars) from the analysis. 

An initial concern when correlating the line shape, 
which is determined from the Civ FWHM, with the 
BH mass residuals, which are also a function of the C IV 
FWHM, is that the correlation is a consequence of this 
redundant dependence, effectively plotting FWHM vs. 
FWHM. To test this possibility, we also calculated the 
C iv BH mass using ui and compare the resulting mass 
residuals with H/3 to the inverted the line shape param- 
eter, S -1 , in the bottom panel of Figure [5] The H/3 SE 
BH mass is still determined using the FWHM. Now both 
axes depend on uciv i n the same way they depend on 
FWHMc iv in the top panel. If the correlation in the top 
panel was due only to a redundant dependence on line 
width in both coordinate axes, we would expect to see a 
correlation of similar strength in the bottom panel, but 
wc find no correlation, with r rs — 0.094 and P ran = 0.533. 
Furthermore, a partial correlation test between the quan- 
tities correlated in the top panel, with FWHMc rv as a 
nuisance variable, results in a partial correlation coef- 
ficient, r — 0.389, which is nearly equal to that from 
the equivalent test for the bottom panel, with crew as 
the nuisance parameter, i.e., r = 0.381, yet again, we 
don't see a correlation of the crciv-based masses with 
S^ 1 . This not only confirms that the strong correlation 
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Figure 5. Top: Comparison between Civ and H/3 mass resid- 
uals, calculated using FWHM, and the shape of the Civ line, 
S = FWHM/o-, for the All (solid circles), VP06 (open circles), 
and N07 (stars) samples. The solid line shows the best fit correla- 
tion to the combined sample, excluding absorbed N07 objects (gray 
stars). The Spearman rank order coefficient, r B , the number of ob- 
jects included in the fit, TV, and the probability that the correlation 
is found by chance, P r an, are given in the top right corner. Error 
bars are not included on the mass residuals because uncertainties 
could not be calculated similarly for all samples (see Section 
Bottom: Same as top panel except the C IV masses were calculated 
using o~ i and the resulting residuals with the H/3 masses are com- 
pared to S 1-1 so that both plots have the same dependence on the 
C IV width characterization used in the mass estimates. Note: r B 
is insensitive to the choice to use S or S^ 1 as the independent 
variable. 

of FWHM-based masses with shape is real (i.e., not due 
to a redundant dependence of BH mass and shape on line 
width) , but also that C iv masses based on 07 more reli- 
ably reproduce, in general, the H/3 mass estimates. This 
latter point is not altogether surprising, since 01 depends 
less on the structure of the core of the line profile, where 
the non- variable C iv emission component most strongly 
affects the shape. 

The strength of the correlation between the FWHM- 
based mass residuals and Civ line shape suggests that 
this shape trend could explain a significant fraction of 
the scatter between SE C iv and H/3 masses observed in 
the literature, as the studies discussed in the Section Q] 
generally used FWHM-based C IV masses. We can fit the 
correlation of the mass differences with shape to derive 
a corrected FWHM-based C iv SE mass equation of the 
form 



log Mqxv (FWHM) = log M°"s (FWHM) + 0.219 

V fa iv / 

where the uncertainties in the slope and the intercept 
are given in Table [5] The top panel of Figure [5] shows 
the original FWHM-based C iv-to-H/3 mass residuals as a 
function of H/3 SE BH mass. The scatter in the combined 
sample about the 1:1 relation of equal Civ- to- H/3 mass 
is 0.39 dex, and there is a mean offset of —0.23 dex. The 
bottom panel shows the same residuals after correcting 
for the Civ line shape bias using Equation (1). There 
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Table 2 

Fit Parameters to Correlation between BH Mass 
Residuals and C IV Line Shape 



Data Set 


Slope 


Intercept 


Scatter 


All 


3.65 ± 0.79 


-0.336 ± 0.101 


0.20 dex 


N07 


1.49 ± 0.67 


-0.172 ± 0.143 


0.16 dex 


VP06 


2.01 ± 0.64 


-0.046 ± 0.062 


0.26 dex 


Combined 


1.63 ± 0.38 


-0.219 ± 0.051 


0.26 dex 
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Figure 6. C lv-to-H/3 SE BH mass residuals from the VP06, All, 
and N07 samples as a function of H/3 mass. Symbol types are the 
same as Figure [5] The top panel shows residuals calculated with 
the original C IV masses. The bottom panel shows the residuals 
after correcting the C IV masses using Equation ( 1 ) . The mean 
scatter about and offset from (dashed line) the 1:1 relation (dotted 
line) for the combined sample is given in the top corner of each 
panel (The corrected scatter and offset exclude the N07 absorbed 
objects.). 

is no longer a mean offset and the scatter is now only 
0.26 dex. We have excluded the absorbed N07 targets 
(gray stars) from the fit and latter scatter determination 
because the C IV mass of these objects is biased by factors 
other than the intrinsic ClV line shape (e.g., absorption; 
see also Denney et al. in prep.). 

Despite the decreased scatter, there are still outliers. 
Many of these are the previously discussed N07 objects. 
The largest All sample outlier is QO 142 -100, but this 
object was classified as a group II object by All, due to 
an unreliable H/3 line width, so the problem may have 
nothing to do with C iv. The two largest VP06 outliers 
(one above and one below the 1:1 relation) are the two 
radio loud AGNs in the RM sample, 3C 120 and 3C 390.3. 
This may simply be a coincidence or may indicate that 
radio loud AGNs follow a different C iv shape relation 
than radio quiet AGNs. For example, iRichards et al.l 
(|2011[) found very distinct differences between the C IV 
properties of the radio-loud and radio-quiet populations. 

5. DISCUSSION 

5.1. Other Line Width and Shape Characterizations 

The FWHM is the most widely-used line- width charac- 
terization in the literature, and it has the advantages of 
being easy to measure and relatively insensitive to blend- 



ing in the line wings. On the other hand, measurements 
of 07 are robust for a wide range of line profiles and have 
been found to be a less biased line characterization in 
some cases (see, e.g.. [Peterson et al.ll2004l iDennev et al.1 
120091 : iRafiee & Hallll2011| ). The problems with measur- 
ing 07 include difficulties in accurately defining the line 
wings, particularly in low S/N spectra, and in the pres- 
ence of other blended emission-line features. Despite 
these concerns, Figures Q] and [5j respectively, show that 
characterizing the C iv width with ai leads to (1) veloc- 
ities that are more consistent with the 1:1 relation be- 
tween the mean and rms spectrum in the reverberation- 
mapped sample, and (2) masses that show less bias with 
the line shape, S. In fact, when we calculate the C iv 
mass using 07 (i.e., Equation 8 of VP06), the scatter 
in the mass residuals about the 1:1 relation of 0.28 dex 
is only slightly larger than the scatter in the corrected 
FWHM-based masses (0.26 dex) shown in Figure El al- 
beit with an additional systematic offset of —0.17 dex in- 
dicating a zero-point calibration difference between the 
full sample and VP06. This offset could be due to the 
C iv shape bias or other systematic uncertainties in the 
exact prescription for calc ulating 07, such as blending 
with the red shelf (see, e.g. iFine et aLll2010b lAssef et all 
[20T1 . In either case, if we adjust the zero-point of 
the VP06 relation to eliminate the offset, the remain- 
ing scatter about unity is only 0.26 dex, equivalent to 
the FWHM-corrected Civ masses. 

Because our shape correction to the C IV mass is a com- 
bination of only the FWHM and 07 , expanding Equation 
(1) shows that we have effectively fit a virial relation of 
the form M oc FWHM 1 of , where y = 2 - x in gen- 
eral, and here, y = 1.63 and x = 0.37. The preferen- 
tially higher weight of the line dispersion in this expan- 
sion explains the similarity of the scatter between the 
shape-corrected masses and cr/-based masses. We also 
determined the Civ shape correction strictly from the 
kurtosis of the C iv profile. Results are not shown be- 
cause they were only based on the N07 and All samples 
for which we had access to the data, but are consistent 
with using S — FWHM/07. This demonstrates that uti- 
lizing either (1) a combination of line width plus shape, 
or (2) a more complex (i.e., than FWHM) characteriza- 
tion of the C iv line profile is necessary to properly cali- 
brate C iv and H/3 masses due to t he contamination of the 
non- variable component of C iv. iBaskin fc Laorl (|2005l ) 
also found that the scatter between C iv and H/3 masses 
could be reduced by fitting a C iv mass correction based 
on additional C iv profile information, albeit with more 
parameters, to wh at we have done here . Additionally, 
iWang et all (|2009[ ) and lPark etai] (|201lD have suggested 
relaxing a strict virial requirement (M oc V 2 ) altogether 
as a means to reduce scatter in the calibration of SE mass 
scales due to systematic differences between various SE 
line profiles and the rms H/3 profil e to which these mass 
scales are ultimately calibrated. IRafiee fc Halll ()2011[ ) 
found that this relaxation could also explain the "sub- 
Eddington boundary" for FWHM-based BH masses, but 
then again, so could using cr/-based masses with the typ- 
ical virial requirement. 

5.2. Sources of Remaining Scatter 

While the ClV shape correction defined by Equation 
(1) significantly reduces (by 0.13 dex) the scatter in the 



C iv-to-H/3 mass residuals, a scatter of ~0.3 dex still re- 
mains. We briefly investigated three potential sources 
for this remaining scatter: 
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Data quality: Several absorbed N07 objects (gray 
stars) became even larger outliers after the C iv 
shape correction. This suggests that the measured 
line width in these objects is not the intrinsic line 
width, so the correction to these sources was mean- 
ingless. This is a serious concern for SE masses 
measured from survey quality spect ra, which are 
typically of relatively low S/N (see IDennev et al.l 
2009; [Asset et al.ll2011l and Denney et al., in prep, 
for additional discussions and investigations of this 
concern) . 

Data Inhomogeneity: As there is no univer- 
sal prescription for AGN spectral decomposition 
and line width measurements, inhomogeneity in the 
data analysis and computational methods used to 
derive SE masses leads to measurable systematics 
when trying to compare masses between samples 
and emission lines. This is demonstrated by the 
comparison of our FWHM mea surements of the 
N07 sample in Figure B) (see also lAssef et a.1.11201 it 
Vcsterga ard et al.l 1201 lL for additional sample to 
sample comparisons). There are also no universal 
mass scaling relationships for each emission line, 
and differences in the mass scale also lead to differ- 
ences in ma sses derived by different studies (see di s- 
cussions by IMcGill et alBool IShen fc Liuil2012l) . 

H/3 Systematics: Literature comparisons be- 
tween SE C iv and H/3 masses typically blame C iv 
masses for any observed inconsistencies, since H/3 
is arguably the most well-characterized and best- 
studied line in terms of direct RM-based masses 
for comparison. However, H/3 it is not without sys- 
tematics of its own, particularly contamination by 
NLR emission and host galaxy starlight, which are 
relatively unimportant for C iv. The systematics 



luminosity can be significant ( 


Dennev et al.l 2009; 


IBentz et al.H2009al IPark et al.l 
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We tested if the scatter in the C iv-to-H/3 masses could 
be further reduced by addressing these additional sources 
of scatter. First, we removed (1) single VP06 epochs 
with S/N < 10 pixel -1 in the continuum, (2) the seven 
objects from the N07 sample that show evidence for ab- 
sorption (we could not cull S/N < 10 pixel -1 data from 
this sample without eliminating it completely), and (3) 
the three All sources defined as Group II objects due 
to un reliable H/3 widths and other issues fsee lAssef et al.l 
[20TTI for details). 

Next, we removed host starlight from the VP06 sample 
by replacing the SE H/3 masses with the RM H/3 masses. 
Ideally, host starlight should be removed from all sam- 
ples, but the required host fluxes are not available. In 
any case, the VP06 sample is the most susceptible to this 
bias because of its relatively lower mean AGN luminos- 
ity. This choice removes some of the consistency between 
the treatment of this sample compared to the All and 
N07 samples, but these differences are smaller than the 
accuracy gained by this correction. 
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Figure 7. Same as Figure [6] except the top panel shows the sub- 
set of data kept from each sample after removing poor quality 
data, removing the All "group II" objects, and using the RM H/3 
masses instead of SE masses for the VP06 sample. The bottom 
panel shows this same "high-quality" sample after applying C IV 
mass corrections similar to Equation (1), but determined for each 
individual sample. Mean sample offsets are listed but not marked 
because they are smaller than the typical measurement uncertainty 
on the mass residuals, ~ 0.1 — 0.2 dex for the VP06 sample. 

Finally, fully removing the data analysis inhomogene- 
ity is not possible. However, we fit and defined a shape- 
based C iv mass correction independently for each sam- 
ple, since the sample-to-sample inhomogeneities in the 
spectral fitting and line width measurement techniques 
have likely led to a sub-optimal C iv shape correlation 
fit based on the combined sample. We list the individ- 
ual sample fit parameters and uncertainties in Table [5J 
Figure [7] shows the results of these modifications to our 
analysis. The top panel again shows the ratio of the C IV- 
to-H/3 masses of this "high-quality" sample before cor- 
recting C IV masses for the shape bias. Here the scatter 
around the 1:1 relation is now 0.33 dex, and the system- 
atic offset is small, showing that simply removing poor 
quality data and host galaxy starlight alone can reduce 
the scatter. The bottom panel shows the residuals after 
correcting the C iv masses using the individual shape cor- 
rections determined for each homogeneous sample. Here 
the scatter is only 0.22 dex, demonstrating the Civ and 
H/3 masses to be in very good agreement. The remaining 
scatter is less than the previously estimated intrinsic sys- 
tematic uncertainties of ~ 0.3 dex inherent in SE masses 
(VP06) and is comparable to the un certainty in the cal- 
ibration of the R — L relation ship ([Bentz et al.l [2009a, 
though see also iPetcrson 20Kj for evidence that the in- 
trinsic R — L scatter is somewhat lower) . 

5.3. The BH Mass - Luminosity Color Correlation 

In their comparison of C IV and H/3-based masses in a 
sample of lensed quasars, All discovered a strong cor- 
relation between the ratio of C iv-to-H/3 masses and the 
ratio of UV-to-optical luminosity. Such a correlation is 
naively expected, since the ratio of these masses is a func- 
tion of the ratio of UV-to-optical luminosity. However, 
the slope of the correlation measured by All was in ex- 
cess of that expected, if the only source of the correlation 
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Table 3 

Spearman Statistics for Correlations with 
UV-to-Optical Luminosity Ratio 



Correlated 
Quantity 


Data Set 


r s 


Pran 


N 


S 


Combined 


0.231 


0.152 


40 


S 


All 


0.555 


0.077 


11 
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logCLO-f/iO^) 


Combined 


0.207 


0.201 


40 


log(M CI v/M H) 9)- 










log(L°f7L° p ?) 


All 


0.527 


0.096 
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was the dependence of Mbh on L (Mbh oa L 1 / 2 ). Inter- 
estingly, if we combine the blue shift-equivalent width re- 
lation of iRichards et all (120111) w ith the blueshift-color 
relation of iGallagher et al.l (|2005| ). we expect a shape- 
color correlation in the sense that high (low) S C IV pro- 
files are seen in bluer (redder) AGNs. This means that 
our shape correction operates in the same sense as the 
All color correction. 

Figure H] shows this shape-color correlation. To be con- 
sistent, we used only SE, non-host-corrected AL5100 for 
all objects, which necessarily decreased the VP06 sam- 
ple by six objects. We compare the statistics from this 
shape-color correlation with the All mass-color cor- 
relation, but we first subtracted the expected depen- 
dence of the All correlated quantities on the ratio of 
the luminosities by removing the the luminosity term, 
(ALi35o)°- 53 /(AL 5 ioo) 52 , from the All mass residuals. 
This leaves only the excess correlation observed by All, 
and is a more consistent comparison with the shape-color 
correlation, in which there is no redundant dependence 
on luminosity in both coordinate axes. The Spearman 
statistics listed in Table [3] show that both correlations are 
weak, but the correlation of S with color is the stronger 
of the two. Since All did not use this combined sam- 
ple to discover or define this correlation, we also make 
a comparison using the All sample alone. The fit to 
the All sample is denoted by the solid line in Figure [51 
and the Spearman statistics are listed the top corner as 
well as in Table |3l The statistics show that these cor- 
relations are still not particularly strong, but are more 
significant than for the combined sample (as also noted 
by All). The correlation of S with color is, again, the 
more statistically significant of the two. 

We also measured the mass-color correlation slope of 
the All Group I objects after correcting these FWHM- 
based ClV masses with Equation (1) above. All origi- 
nally measured this slope to be 0.89 ±0.25. We measure 
a new slope of 0.71 ±0.17, which is marginally consistent 
with the expected slope of ~0.53 from the mass scaling 
relationship dependence on L. This suggests that the 
shape bias in the Civ masses is a likely driver for the 
All mass-color correlation and can at least partially ex- 
plain the steep slope observed by All. Additional work 
is underway using a larger sample to either identify or 
rule out additional effects, e.g., non- universal SEDs, that 
could also be contributing to the excess slope observed 
by All. 
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5.4. Origin of the Non-variable C iv Component 
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Figure 8. Comparison between Civ line shape, S = FWHM/cr;, 
and the ratio of the UV-to-optical luminosities, ^1350/^5100 = 
AL A (1350A)/AL A (5100A). Symbol types arc the same as in Figure 
[5] The solid line is the best fit to the All sample alone. The 
statistics in the top left corner are for the All sample only; the 
combined sample has r B = 0.231 and P ran = 0.152 (see Table|3}. 

Our analyses in §S}3] and [4] have shown (1) that the SE 
C iv profile is a composite of both a variable (reverber- 
ating) component and non-variable component, and (2) 
that object-to-object differences in the characteristics of 
the non-variable component are a primary cause for the 
large scatter in FWHM-based SE C iv BH masses as com- 
pared to H/3 masses. Here we discuss possible origins for 
these differences. 

5.4.1. The "Traditional" Narrow Line Region 

The first, obvious, possibility is to associate the low 
velocity component with the "traditional" , low-density 
NLR, i.e., the region responsible for emitting narrow for- 
bidden lines, such as [O III] AA4959, 5007. A NLR com- 
ponent is observed in the H/3 emission line and has some- 
times been as sumed to be a component of the C IV profile 
as well (e.g. IBaskin fc Laor|[2005l: ISulentic et alf l2007l: 
iGreene et all 120101 : iShen k. Liul 12012ft . However, we ar- 
gue that this cannot be the origin of the non-variable 
component of C IV for the following reasons: 

1. While the C+3 ion is photo-ionized, the C iv A1549 
emission line transition is excited collisionally. 
Photoionization models suggest that Civ A 1549 
emission from the l ow-density NLR region mus t 
therefore be weak (jFerland fc Osterbrockl [1 986) . 
In contrast, the non- variable component seen in 
the mean C iv profiles of Figure [3] is often quite 
strong. Admittedly, Civ e mission is sti ll seen 
in Seyfert 2 spe ctra (e.g., iCollins et al.l 120051 : 
iZheng et al.l 12008) , also suggesting a possible ori- 
gin in the NLR, but studies have additionally 
shown that NLR emission line strengths and ra- 
tios are different betw een Seyfert l's and Seyfert 
2's dZhang et al.ll2008|) and may be dependent on 
inclination (|Fine et al.l 1201 If) . 

2. Not all Civ SE profiles show an obvious nar- 
row, core component. The presence and strength 
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of this component is also anticorr elated with 
blueshift and the C IV equivalent width (iWills et all 
1991 IRlchards et al.ll200l ILeighlv fc Moorell2004t 
Richards et all 1201 ll ). It is much more probable 
that this reflects differences in the BLR kinematics 
or geometry between objects rather than correla- 
tions between NLR and BLR kinematics. 

A large per centage of Civ SE p rofile peaks are 
blueshifted (|Richards et all 120111 ). some by thou- 
sands of km s _1 . While NLR emission can be 
blueshifted from the systemic redshift, it would 
not be expected to show these large of blueshifts. 
If present, NLR emission should appear at or at 
least near the systemic redshift, leaving it relatively 
redshifted with respect t o the highly blueshifted 
Civ profile (Wilkes 1984). One would then expect 
potential double-peaked profiles in these objects, 
rather than this component simply being absent. 
Such a feature is not observed here, in large sam- 
ples, or even in com posite AGN spectra (see, e.g., 
IRlchards et aHl20ll . 

We subtracted the scaled rms spectrum from the 
mean spectrum of all objects shown in Figure [3J 
(except Fairall 9 because of the poor quality rms 
spectrum) and measured the FWHM of the resid- 
ual non-variable components using the same meth- 
ods as in Section |31 We also measured the FWHM 
of the O in A5007 line from this sample, subtract- 
ing a linear, local continuum, from available optical 
data. Figure[S]shows that the velocity widths of the 
non-variable narrow cores are larger (often signif- 
icantly so) than those measured in the unblended 
NLR OiiiA5007 line. Combining these relatively 
broad core widths with the core strength and addi- 
tional blue excess emission in most of this sample, 
it is unlikely that this component is, at least solely, 
from the traditional, low-density, extended NLR 
emission. 



5.4.2. A BLR Disk-wind 

The presence of a BLR wind has been used to ex- 
plain a wide variety of observed AGN phenomena. This 
includes, for instance, the of ten blueshifted Civ pro- 
file (see IRichards et all 120111 a nd references t h erein) , 
absorption line features (e.g., IHamann" 
and A GN outflows and feed back (e.g., iReeves et all 
2009]). IRichards et all (|20ll argue that the sources 
with the largest blueshifted Civ profiles are the 
"wind-dominated" sources, where this wind component 
causes a blueshifted peak and line asymmetries (cf. 
ILeiehlv fc Mooril200l in the SE C iv profile. These ob- 
servations are consistent with the non-variable blue ex- 
cess emission we observe when comparing the RM mean 
and rms spectra. However, we also observe a strong low- 
velocity, non-variable core in some objects. So in terms 
of the SE BH mass estimates alone, if this total non- 
variable component (core + blue excess) is a disk-wind, 
it is the peaky, low S- value objects with the strongest 
SE core components that would have the largest wind 
"contamination" or bias in their FWHM and therefore 
FWHM-based SE mass measurements. 




_i 1 1 1 1 1 1 1 1 1 1 1 1 1 1 

500 1000 1500 2000 

[0III]\5007 FWHM (km S"') 

Figure 9. Comparison between the FWHM of the C IV non- 
variable component and the FWHM of the [Oiii] A5007 line for 
the RM sample. Individual objects are labeled, and the measure- 
ments reflect all objects shown in Figure [3] except Fairall 9, whose 
rms spectrum was too noisy to accurately isolate the non-variable 
Civ component. The solid line is the 1:1 relation in FWHM, and 
the two Civ non-variable FWHM measurements for NGC 4151 
connected by the dotted line represent the difference in width if 
the narrow absorption near 1545A is linearly interpolated across 
before measuring the FWHM (upper value) or not (lower value). 

The BLR disk-wind has been suggested to be launched 
from the inner, high-velocity, high-io nization regions of 
the BLR or outer accr etion disk (cf. iMurrav fe Chiang] 
119971: lElvisI 120001 12012D . so the additional non- variable, 
low-velocity core component observed here does not read- 
ily fit into this model unless there is an orientation depen- 
dence of the wind component relative to the line of sight 
(LOS). The RM mean and rms spectra in Figure [3J as 
well as C iv profiles in general show that the strength of 
the blueward asymmetric emission appears to be roughly 
anti-correlated with the strength of the low- velocity core 
(e. g., the blueshift-equi valent width relation discussed 
bv lRichards"et~aLl l20ll . Therefore, we argue that if 
this non-variable component does originate in an out- 
flowing BLR wind, there should be an orientation depen- 
dence, where the change in C IV profile from 'peaky' to 
'boxy', or low S to high S, is due to varying levels of 
LOS wind contamination to the broad emission line pro- 
file as a function of orientation. Furthermore, because 
this component does not reverberate, such a wind would 
necessaril y be optically thin to the ionizing continuum 
(see, e.g., iKorista fc Go ad 2004, who discuss the possi- 
bility that low optical depth can explain the non- variable 
wings of the H/3 profile). 

For a disk+wind BLR (cf. Figure 1 of 
IMurrav fc Chiang] 119971) to support our observa- 
tions of the variable+non- variable C iv profile and put 
these observations in context to other C iv observations, 
an optically-thin wind that is more polar than equatorial 
is required. In this case CivA1549 can still be emitted 
in the wind but does not reverberate. We can then make 
a simple generalization for a composite "disk" (variable 
component) plus "w ind" (non-variabl e component) 
Civ profile (see also iWang et ail 1201 ID . We would 
then expect the following orientation-dependence of the 
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profile of each component: 

• The disk: As the inclination decreases (goes from 
edge-on to face-on), the disk component experi- 
ences a general narrowing in velocity width due 
to the change in inclination. There is no other a 
priori requirement for the shape of the line; pho- 
toionization and BLR geometry model-fits to the 
Balmer line profiles have shown that a wide variety 
of both single- and double-peaked profiles can exist 
and still arise from a virializcd disk-like distribution 
of gas with varying o ptical depths (see review by 
lEracleous et al.ll2009l and references therein). 

• The wind: At high (closer to edge-on) inclina- 
tions, the wind component is narrow and centered 
near the C iv systemic velocity, as the direction of 
outflow is largely perpendicular to the LOS. As the 
inclination decreases, more of the wind is observed 
to be outflowing along the LOS. Consequently, the 
wind emission is distributed across a larger, pref- 
erentially blue-shifted, range of velocities, so the 
peak also blueshifts and decreases in flux. 

5.4.3. The Intermediate Line Region 

Another previously suggested explanation for the 
peaky SE core component in some C IV lines is that this 
emission arises in an "intermediate line region" (ILR; 
iWills et al.l 119931: iBrotherton et all 119941) . Such a re- 
gion is explained as a higher-velocity, higher-density in- 
ner (i.e., c l oser to the BLR) extension of the NLR. 
IWills et all (|1993H alternately suggest a different, simple 
model in which the core C iv emission could be coming 
from a bipolar outflow, but this may not be mutually 
exclusive to the ILR/inner NLR. More recent spatially- 
resolved studies of the kinematics and ph ysical conditions 
in the inner NLRsof local Sevferts (e.gjDas et al.ll2005t 
Kraemer et alJl2009l iCrenshaw et alJl2010l : lFischer et al.1 
20101120111 1 show a link between p ossible ILR C IV emis- 
sion and biconical outflows (cf. Crens haw fc Kraemer! 
2007). While we mentioned above that densities in 
the "typical", i.e., extended, NLR are not high enough 
to result in significant emission from high-ionization 
species, spatially resolved STIS spectra have co-located 
high-ionization-line emission with knots of [O III] A500 7 
emission in the inner NLR (e.g. ICollins et all 12005). 
Such e mission must be enhanced by collisional pro- 
cesses (|Kriss et al.lll992T) . such as shock heating or mi- 
croturbulance, mechanisms likely to be present if thi s 
emission is arising in a n outflow dKraemer et al.1 12007( 1 . 
iNelson et all ((2000( 1 and lKraemer et al.l ((2000( 1 show that 
with HST observations of nearby Seyferts, spatially- 
resolved Civ A 1549 emission does arise outside the un- 
resolved nucleus, i.e., from the ILR/inner NLR. How- 
ever, the C iv emission flux drops significantly faster 
than [O in] AA4959, 5007 emission as a function of ra- 
dius. The spatially resolved inner knots of [Olli] A5007 
emission also show relatively high, ~ 1000 km s -1 , and 
even blueward a s ymmetric velocities ((Nelson et al.ll2000t 
iDas et all (20051 : iFischer et all 1201 11 1, both reminiscent 
of the non-variable C iv component we see here. Un- 
fortunately, and important for potentially characterizing 
the non- variable C IV profile, because the [O in] emission 
extends to much larger radii and lower velocities, the 



line profile of the "integrated" [O in] emission in typical 
ground-based spectra of both near and distant objects is 
not a good template match for the observed non- variable 
C IV core emission. 

In this type of scenario, the non-variable component 
of the Civ line, i.e., both the peaky, low-velocity Civ 
core emission and blue excess, could be attributed to 
the ILR emission and part of a biconical outflow. LOS 
emission-line profiles from such a distribution of gas 
would therefore exhibit velocity widths larger than the 
typical, extended NLR region that is probed by in- 
tegrated forbidden-line emission profiles, but narrower 
than typical BLR emission. Additionally, if the gas is 
outflowing, blueshifts are to be expected, and poten- 
tially more prominent for the higher gas velocities, since 
this emission likely arises at closer radii and may there- 
fore by more highly collimated by the bicone. These 
are all expectations consistent with the observations of 
the C IV emission line profile seen here and in previous 
studies. Furthermore, our orientation-dependent model 
of the non- variable C iv emission introduced above fits 
equally well with an origin in a biconical ILR outflow as 
in the BLR disk-wind, since the LOS component from the 
inner region of such an ILR outflow should change as a 
function of orientation. Such a distribution of gas would 
also not be expected to reverberate. Because the spa- 
tial scales are large relative to the time scale of variabil- 
ity, any reverberation signal would be spatially damped. 
Therefore, the non-variable nature of such an emission 
component naturally holds for this origin, at least on re- 
verberation time scales. 

5.5. Evidence for Orientation- Dependence 

Both origins for the non- variable C iv component that 
we suggest above require an orientation dependence to 
the outflow component of C iv emission. We first looked 
for indications of an orientation dependence by compar- 
ing our orientation expectations with the Civ profile 
trends observed in the RM sample. To do this, how- 
ever, we must isolate the non- variable component profile 
to search for signs of orientation dependence. This was 
done through a spectral decomposition of the mean spec- 
trum C iv profiles shown in Figure [3] We assume that 
the mean profile is a superposition of the BLR 'disk' 
(variable component) and outflow (non-variable compo- 
nent). The rms profile shows variable emission only, so 
we generalize this emission to arise solely in the reverber- 
ating BLR 'disk'. We decompose the mean profiles into 
variable and non- variable components by first fitting the 
scaled rms profile (solid gray curves in Fig.[3|) with a sixth 
order Gauss-Hermite polynomial as a means to derive a 
smooth profile for the disk component. We subtract this 
component from the mean spectrum, and attribute the 
residual non- variable flux to the outflow component. We 
fit the residual emission again with a GH polynomial af- 
ter interpolating over any narrow absorption features. 

We highlight this decomposition for the Civ profiles 
of NGC 4151, whose mean spectrum has a relatively low 
iS-value (0.61) with a strong non-variable core, and NGC 
7469, which has higher S (0.92) with a notably weaker 
non-variable core but with additional, non-variable blue 
excess. The results of the decomposition for these two 
sources are shown in Figure 1101 where the original mean 
and rms spectrum and individual and composite (disk 
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Figure 10. Predicted emission line profiles expected if observing 
an object with a disk-wind BLR at higher, more edge-on inclina- 
tions, (top) and lower inclinations (bottom). For each panel, the 
higher flux black curves show the original mean spectrum (also 
shown in Fig. [3]l used to represent fiducial C IV profiles for each 
orientation. The red curve represents the variable model compo- 
nent fit to the scaled rms spectrum (also in black), the blue curve 
represents the non- variable model component, and the gray curve 
is the composite model including both components. Finally, the 
dashed curve is the AGN continuum that was subtracted before 
fitting each emission-line component. 

plus outflow) profile fits arc shown for NGC 4151 and 
NGC 7469 in the top and bottom panels, respectively. 

From this exercise, we observe for NGC 4151 that (1) 
the rms profile (our disk component) is broad and double- 
peaked, possibly dem onstrating the emission exp ected 
from an inclined disk (|Eracleous fc Halpernl 17994) . and 
(2) the relatively narrow velocity distribution of non- 
variable residuals would also be expected from an out- 
flow directed largely perpendicular to the LOS. While 
the model-fit to the rms spectrum over-predicts the flux 
in the core, this acts only to somewhat underestimate the 
core strength but does not significantly change the pro- 
file of the residual non- variable component. As it would 
not be possible to view this broad-line AGN perfectly 
edge-on and still see a broad-line AGN (in the unified 
model), the extended and blue-shifted wing in the non- 
variable component profile is also consistent with the ex- 
pectations for this component arising in an outflow. For 
NGC 7469 we see (1) a less peaky and more blueshifted 
outflow component and (2) a narrower disk component. 
These are both expected characteristics of a more face- 
on orientation. The combined fit conspires to produce a 
smooth but slightly asymmetric profile, characteristic of 
many SE C iv profiles. 

The orientation-dependent model we describe here is 
naively simple. Interestingly, however, our crude decom- 
position using two test cases from the RM sample places 
no a priori expectations for the profiles of the individual 
variable and non-variable components. Yet, the results 
of attributing these components to a BLR disk+outflow 
composite profile are strikingly similar to what is ex- 
pected if we are viewing NGC 4151 at a higher inclina- 



tion (more edge on) and NGC 7469 at a relatively lower 
inclination (more face on). 

There are also numerous discussions in the liter- 
ature of orientation depe ndencies to ob s erved AGN 
properties. For example, iRichards et al.l J2002 : ) sug- 
gested the possibility of an orientation effect in the 
Civ blueshifts, and thus wind component, though in- 
terestingly in t he op posite sense that we argue here. 
IRichards et al.l ()2002f ) advocated that large blue-shifted 
profiles were viewed from more edge-on orientations, 
which follows t he expectations of the l argely equato- 
rial w ind in the IMurrav &: Chiang! (|1997t ) model. How- 
ever, IRichards et al.l ( 20111) argue against such a de- 
pendence based on radio morphology studies. Other 
radio morphology studies also show no or only weak 
orientation-dependenc e of C iv widths ([Jarvis fc McLurel 
120061 ; iFme et al.ll201lD . However, these studies are based 
primarily on using the FWHM, and we have shown here 
that C iv FWHM measurements are biased by the non- 
variable component. In fact, Civ profiles resulting from 
an orientation-dependence of the non-variable compo- 
nent in the sense we propose above (creating preferen- 
tially peaky profiles for high inclination and boxy pro- 
files for low inclination) would act to mask the expected 
changes in the FWHM of the BLR 'disk' component that 
these studies hope to probe: The non-variable compo- 
nent contamination will lead to an artificially large mea- 
surement of the FWHM for low inclinations and an arti- 
ficially small measurement of the FWHM for high incli- 
nations. 

Another place to search for possible orientation- 
dependent quasar phenomena is in the characteristics of 
broad absorption line quasars (BALQSOs), which have 
been inferred to be viewed preferen tially more edge-on 
(e.g.. IMurrav et al.l[l995t lElvisl[200 01. This is partly be- 
cause their conti nua are relative l y mo re reddened than 
non-BALQSOs (|Reichard et al.l 120031 and references 
therein). However, other observations and models have 
shown that (1 ) polar outflows a re als o capable of pro- 
ducing B ALs (IBrotherton et al.l 120061: iZhou et al l 120061 : 
iGhosh fc Punslvl l2007t iBorguet fc Hutsemekersl I2010T ). 
(2) individual fits to BALs cannot c onstrain the geometr y 
and orientation of the BAL region (jHamann et al.f l993). 
and (3) radio BAL and nonBAL quasars are nearly indis- 
tinguishable across a wide range of observed and physi- 
cal properties, including sp ectral shape, spec tral index, 
and polarization properties (Br uni et al.ll2012l ) . This ev- 
idence suggests that BALs can be observed in quasars 
over broad range of orientations. 

Furthermore, discussing the results of some previous 
studies in the context to this simple orientation depen- 
dence of the C IV SE profile leads to some novel interpre- 
tations of the observations (though we do not claim to be 
able to explain every aspect of these observed phenomena 
with only inclination) : 



The Baldwin Effect (jBaldwinl Il977l ). We argue 
that the lowest equivalent width objects are those 
that tend to have high 5* and are the most face-on. 
The continuum photons escaping the central source 
should therefore be the least obscured by other nu- 
clear material. This will lead to a larger observed 
luminosity, on average, compared to higher inclina- 
tion objects, which we argue are observed to have 
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high-equivalent width, low S profiles due to wind 
contamination. 

• NLR Component Remov al (e.g. iBaskin & Laorl 
120051 : iSulentic etaD 120071 ). Some studies argue 
the narrow C i v core originates in the NLR. Yet, 
ISulentic et al.l (|2007l) have pointed out that the 
narrow-line components in C iv are often broader 
and stronger than typical, "Qui A5007 -like" , NLR 
emission. Though ISulentic et al.l (|2007t ) found bet- 
ter correlations between C iv and H/3 after remov- 
ing a narrow core component, they were still left 
with a large scatter between their C iv and H/3 
masses. In the context of our model and other argu- 
ments in Sectio n 15.4. II the large scatter remaining 
in the ISulentic et a l. ( 20071) results is most proba- 
bly due to the uncertainty in fitting and subtracting 
their "narrow" component. Since the orientation- 
dependent non-variable component cannot simply 
be described as a typical NLR emission line, a tra- 
ditional spectral decomposition does not work for 
C IV. 

• "Extincte d" red wing of ClV (see e.g ., composite 
spectra of lRichards et all 120021 I2011 . Pa st phys- 
ical e xplanations for this phenomenon (cf. iGaskelll 
1982) include a combination of radial gas motions 
and obscuration, though Gaskell admits that it is 
questionable whether the dust grains required for 
such obscuration could survive in the BLR. Our 
orientation-dependent disk-outflow model explains 
this phenomenon simply: as the inclination de- 
creases, the observed range of disk velocities nar- 
rows. Conversely, the observed velocities of the 
outflow component are increasing as it becomes di- 
rected more along the LOS, effectively filling in the 
high-velocity portion of the profile previously oc- 
cupied by the disk emission. However, this process 
is preferentially blueward asymmetric because the 
receding component of the outflow is not visible 
through the disk. This leaves the red side of the 
profile devoid of emission from the highest velocity 
gas. 

• Correlation o f blu eshift with luminosity 
()Richards et al.l 120111) . This is expected if 
these are the population of objects see n more 
face-on. Furthermore, Gallaghe r et al.l ([2005) 
demonstrated with a color-blueshift distribution 
that large blueshift AGNs are less likely to have 
red continua, regardless of whether intrinsic or 
dust-reddened. This would similarly be true if we 
are seeing these objects more face-on. 

• Intrin sic colors of BALQSOs (Rcichard et al.1 
2003 ). Despite the observed trend of BALQ- 
SOs being mor e redd ened than nonBAL quasars, 
iReichard et al.l (|2003l ) found that after correcting 
for SMC-like dust reddening, BALQSOs appear 
to be drawn from the same parent p opulation as 
nonBALs (see also iBruni et^l2012h . Upon fur- 
ther detailed inspection of B ALQSO propert i es as a 
function of 'intrinsic' colors. IReichard et al. (2003) 
also discovered the same nonBAL quasar C iv emis- 
sion line trends of larger (smaller) equivalent width 



profiles for 'intrinsically' red (blue) HiBALs (their 
Figure 9). This suggests to us that the obscuration 
typical of BALQSOs may not necessarily be an ori- 
entation dependent effect (e.g., from a torus, which 
would indicate a preferential edge-on orientation). 
Instead, if the C IV line shape trend is an orien- 
tation effect as we suggest, the obscuring regions 
responsible for producing BALs in quasars may in 
fact be orientation independent. The range of ob- 
served C IV profiles in BALQSOs is then additional 
evidence that BALs can be observed over a range 
in orientations. 

Realistically, there are many geometrical, kinematical, 
and physical properties of the BLR that could cause 
varying contributions from the LOS component of a 
BLR disk-wind, and most likely, more than one effect is 
present in va rying degrees for di f ferent samples of AGNs. 
For example. iMurrav fc Chiang] (|1997f ) find that stronger 
core emission from an equatorial wind component can 
occur without requiring an inclination-dependence; in- 
stead, it is the result of simply extending the outer radius 
of the wind-emitting region. We focus only on inclination 
here because it is arguably the simplest to conceptual- 
ize, still suffices to explain many observed trends, and 
fits the observations whether the non-variable emission 
originates in an outflow from the BLR or ILR. Ideally, 
we would like to use direct orientation measurements 
to corroborate our interpretation that the C iv profile 
changes are a result of an orientation-dependence of the 
non-variable emission. Unfortunately, such orientation 
measurements are scarce. 

The mo st promisin g met h ods use radio morphol- 
ogy (e.g. iBrothertonl 1 19961: IVestergaard et all 120001 : 
Uarvis fc McLurd 1200a Fine et all 120111) or NLR, kine- 
matics (see, e.g.. iFischer et all 120101 1201 If ) . As men- 
tioned above, past radio studies may be biased by typi- 
cally only using the F WHM, and these studies should be 
revisited (though see IVestergaard et all 120001 who find 
the strongest C IV orientation dependence using the full 
width at 20% — 30% peak intensity, which would be less 
susceptible to contamination from the non- variable core) . 
NGC 3783 and NGC 4151, two objects in our RM sam- 
ple, have been targeted for NLR kinematical modeling 
to determine orientation. The best-fit models indicate 
an inclination of only 15° away from the LOS for NGC 
3783 (T. Fischer, priv. comm.; Fischer et al. in prep) 
and 45° for NGC 41 51, near the max imum expected for 
a Type 1 spectrum (|Das et al.l 12005). Putting the rela- 
tive strength of the non- variable core observed in each of 
these objects in Figure [3] in context to these independent 
orientation estimates seems to support our expectations 
that we are observing NGC 3783 and NGC 4151 at rel- 
atively low and high inclination, respectively. 

6. SUMMARY AND CONCLUSION 

We have investigated sources for the large observed 
scatter between C iv- and H/3-based single-epoch AGN 
black hole masses. Our primary line of investigation was 
to return to the assumptions on which SE BH masses are 
based by comparing the C iv profiles in the mean and rms 
spectra from reverberation mapping experiments. These 
include the requirement that the velocity width of the 
mean, or SE profile be a reasonable representation of the 
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rms profile width — the velocity profile corresponding 
to emission from reverberating gas. We found that the 
SE or mean C iv emission-line profiles often have a non- 
variable low- velocity core component and/or blue excess 
that does not appear in the rms profiles and therefore 
does not reverberate. The presence of these non- variable 
components in the SE C iv profiles are biasing C iv-based 
BH mass determinations. 

We developed an empirical correction for this bias, ap- 
plicable to FWHM-based C iv masses by using the ob- 
served characteristics of the C iv profile, which we pa- 
rameterize through the ratio, S — FWHM/07, of the 
C iv FWHM to (Ji . After applying this empirical "shape" 
correction to the FWHM-based C IV SE BH masses in a 
sample of ~ 50 low- and high-redshift AGNs with both 
C iv and H/3 mass estimates, we found that the scatter in 
the ClV-to-H/3 mass residuals is reduced from 0.39 dex 
to 0.26 dex for this somewhat heterogeneous sample. 

We then investigated the remaining scatter in the C IV- 
to-H/3 mass ratio by looking at other possible sources. 
We concluded that data inhomogeneity, poor data qual- 
ity, and independent biases in the H/3 masses due to NLR 
emission and host galaxy starlight also contribute to the 
observed scatter. If we minimize these systematic prob- 
lems, we found that we could further reduce the scat- 
ter in the C iv-to-H/3 mass residuals to 0.22 dex. The 
significant additional decrease in this scatter as a result 
of applying individual shape correction factors to each 
smaller, more homogeneous sample suggests that a shape 
correction calibration determined from a highly homoge- 
neous sample of both C IV and H/3-based SE masses may 
further reduce the mass residual scatter compared to that 
seen here. We are currently in the process of collecting 
such a sample of high S/N, simultaneous observations 
of both the C iv and H/3 emission line regions in a large 
sample of AG Ns using the X shooter spectrograph on the 
VLT (see e.g. JHo et al.ll2012l for a comparatively smaller 
sample of such data). 

Finally, we investigated whether the C IV shape trend 
could also be responsible for the mass-color correlation 
described by All. We found only a weak correlation 
between S with this UV-to-optical color term. However, 
this weak shape-color correlation had a similar slope but 
was actually tighter than the All mass-color correlation 
after removing the expected dependence of that corre- 
lation on L, implying that the mass-color correlation is 
simply a manifestation of the C iv shape bias. We further 
supported this interpretation by demonstrating that the 
steep slope observed by All in the mass-color correlation 
largely disappears after applying the shape-correction de- 
fined in Equation (1) to the All FWHM-based masses. 

Our shape parameterization of the C iv line profile and 
correction to the FWHM-based C iv masses was strictly 
empirical in construction, and no t altogether dissimila r 
from some previous studies (e.g... iBaskin fe Laorl 120051 ). 
However, the additional results gleaned from the rever- 
beration database of mean and rms Civ spectra led to 
the discovery that the observed range of SE C iv line 
shapes correlates with the presence of and object-to- 
object differences in a non- variable emission component. 
The physical interpretation we suggest could explain the 
observed differences in this non- variable component is via 
orientation-dependent sight-lines of a (largely) polar out- 
flow. The SE C iv emission-line profile predictions from 



such a model (peaky in more edge-on sources, boxy and 
blue-shifted in more face-on sources) explains the distri- 
bution of observed SE profiles as composites of both the 
reverberating BLR emission and "contamination" from 
the orientation-dependent non-variable outflow compo- 
nent. This contamination biases the velocity width mea- 
sured in the SE spectrum relative to the intrinsic velocity 
dispersion of the reverberating gas, thus biasing SE mass 
estimates. 

The radial origin of such an outflow remains an unan- 
swered question, however. Profile and velocity width 
comparisons of this non-variable component with the 
[O in] A5007 emission line suggest it must be arising from 
closer in than the extended NLR probed by integrated for- 
bidden line emission. However, it is difficult without fur- 
ther observations and advances in modeling to conclude 
whether the outflowing gas producing the non-variable 
C iv emission originates in an inner BLR disk- wind, the 
intermediate line region, aka, high-velocity inner exten- 
sion of the traditional NLR, or even a combination of 
the two. At present, we show no preference as to what 
that origin may be. Certainly an origin in the inter- 
mediate line region produces natural mechanisms to ex- 
plain the lack of variability of this component (spatial 
damping) and the orientation-dependence (collimation 
in a biconical outflow). Plus, spatially resolved bicon- 
ical outflows have been shown to exist on these scales. 
On the other hand, such an outflow must be launched 
from somewhere within the unresolved nucleus. Models 
seem to favor the inner (high-ionization) regions of the 
BLR or outer accretion disk for the ultimate source of 
AGN outflows/ winds, but the orientation-dependence of 
such a wind is, at present, mo del- dependent. 

Future observational studies could potentially con- 
strain the physical origin of this component by studying 
any potential time dependence of the equivalent width 
and velocity profile of this component. For example, 
gross changes in the broad emission-line profile are ob- 
served to exist over dynamical time scal es for the size 
of the BLR (see, e.g.. iSergeev et al] l2007'). So, while we 
have concluded that this outflow component does not re- 
verberate, if the origin of this component is from an inter- 
mediate line region, variability may still be expected on 
much longer time scales than a single reverberation cam- 
paign. If present, a measure of this variability time scale 
could put a constraint on the location of this emitting 
region. On the other hand, if the emission is originat- 
ing in an optically thin wind close to the central source, 
changes in the amplitude of emission may be expected 
with changes in the AGN luminosity state. As the ioniz- 
ing continuum luminosity increases and the ionization 
front moves outward, we may expect an increase the 
equivalent width of this emission because the volume and 
temperature of the ionized gas will also increase. 

Regardless of the physical origin, this non-variable 
component is likely present in some strength and form 
in all SE C iv emission lines, and we have shown it to 
be a large, but correctable, source of systematic bias in 
C iv-based SE BH mass estimates. While we suggest 
the possibility that the profile of this non-variable 
component is orientation dependent, more sophisticated 
models of the BLR and AGN cent r al engines (cf. the 
mode ls of iChiang &: Murray! [T996t iMurrav fc Chiang! 
[19971 iProgaetal.l 120001 : iProga fc Kallmanl 120041: 
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lOhsuga fe Mineshigd 120111: iProga et al.l 12012ft taking 
into account not only orientation, but also photoioniza- 
tion physics, kinematics, geometry, and even radiation 
pressure are needed. Such models will help us to better 
understand the full complexity and origin of the C iv 
emission (both reverberating and non-reverberating) as 
well as other observed AGN properties. Quasar evolution 
and/or SED and L/L ec jd changes with luminosity could 
also affect the relationship between the variable and 
non-variable C iv emission components. Such changes 
could be independent of orientation and/or lead to 
different orientation dependencies as a function of these 
parameters, if, for instance, these physical properties 
effect the bicone opening angle or the driving mechanism 
for AGN outflows. If such models and/or additional 
observational tests can corroborate inclination as the 
source or at least a large contributor to the differences 
in the observed C IV profiles and for either all or even 
particular samples of AGNs, having an inclination 
indicator as simple as the observed SE C IV profile could 
be a very powerful tool for better understanding the 
physics of AGNs and reducing systematics in both SE 
and RM masses. 

As an interesting side note, if a source of this non- 
variable C iv emission is a BLR disk-wind, we may 
also consider the possibility that the non-variable, high- 
velocity wings of H/3 that disappear in rms spectra and 
have also been potentially attributed t o emission from 
optic ally-thin gas (see discussion by iKorista fc Goadl 
2004), may also be from this same BLR wind compo- 
nent. H/3 is emitted from larger BLR radii than C iv, so 
it may only be the very highest velocities of H/3-emitting 
gas that are affected by this wind component, particu- 
larly as it has been theorized to be launched from the 
inner BLR. 

In the end, however, the rms spectra derived from re- 
verberation mapping campaigns are the only guaranteed 
means by which to directly probe the uncontaminated, 
variable profile most relevant for determining accurate 
black hole masses and to decompose the non-variable 
component by which we can study its origin. There- 
fore, more reverberation mapping experiments targeting 
objects across the full observed distribution of SE C iv 
profiles as well as multiple epochs for individual objects 
are needed (1) to expand the sample of C iv mean and 
rms spectra, (2) to form a more robust C iv shape correc- 
tion, and (3) to better understand the extent and physi- 
cal origin of this non- variable component and the bias it 
imposes on C iv-based BH mass estimates. Furthermore, 
these programs must also probe larger AGN luminosities 
and higher redshifts so that results can more directly 
be applied to the typical quasar populations observed in 
surveys. Currently, our only option is to extend what 
has been learned for the relatively, low-luminosity, local 
Seyfert sample presented here. However, it is imperative 
for our understanding of galaxy evolution to directly test 
what systematics, if any, we are introducing by extending 
results from these local investigations to AGNs at higher 
redshifts, luminosities, and even in potentially different 
evolutionary states. 
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